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METHODE POUR MODELISER EN 2D OU 3D UN MILIEU HETEROGENE TEL QUE LE SOUS-SOL DECRIT 
PAR PLUSIEURS PARAMETRES PHYSIQUES. 

(5^ - Methode pour obtenir un module representatif tf un 
milieu het6rogene en 2D 

OU 3D, decrit par un ou plusieurs parametres physiques 
(tels que les impedances du sous-sol relatlvement a des on- 
des de type P ou S, sa density, etc.), a partir de donnees 
diverses: enregistrements sismiques, diagraphiques, geo- 
logiques, etc. 

- On part d'un modele a priori a un ou plusieurs parame- 
tres forme ^ partir de valeurs ponctuelles connues ou esti- 
m^es d'un ou plusieurs parametres physiques et on 
ameliore le module initial en minimisant une fonction cout 
globale representative d'une part des ecarts entre des enre- 
gistrements synthdtiques tradulsant la reponse du modele 
courant avec les enregistrements sismiques effectu^s, et 
d'autre part les hearts entre le module courant et le module 
a priori. Ces ecarts sont mesures au sens de normes dedui- 
tes des informations a priori sur les Incertitudes relatives 
aux donn6es et au modele, ces Informations pouvant varler 
tf un point a un autre du module. 

- Applications a la caract6rlsatlon de gisements d'hydro- 
carbures par exemple. 
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L' invention conceme une methode pour obtenir un modele 

V 

representatif en 2D ou 3D d'un milieu heterogene tel que le sous-sol, 
decrit par plusieurs parametres ou grandeurs physiques, a partir de 
5 donnees obtenues par des mesures indirectes obtenues depuis Text^rieur 
de celui-ci, et d'autres donnees, notanunent un ensemble de donnees 
ponctuelles mesurees in situ, toutes ces donnees constituant T information 
connue a priori sur le milieu. 

La methode selon 1' invention trouve des applications dans des 
10 domaines tres divers, selon le type d^ondes (ondes elastiques, ondes 
electromagnetiques, etc.) emises pour obtenir ces mesures indirectes. 

Dans le cas d'un milieu tel que le sous-sol, les parametres physiques 
descriptifs sont par exemple Timpedance qu'il presente relativement k 
des ondes P ou a des ondes S ou encore sa densite. Les donnees directes 
15 et indirectes dont on dispose sont par exemple des mesures 
diagraphiques, des donnees d' exploration sismique, et d'autres donnees 
que Ton possMe suite a des etudes geologiques prealables du milieu. 

Les ondes emises dans le milieu par une source sismique se 
propagent dans le sous-sol et se reflechissent sur les discontinuites du 
20 milieu. EUes sont enregistrees par des capteurs couples avec la formation 
souterraine et collect^es par un dispositif d' acquisition. Les donn6es 
sismiques que Ton exploitent dans le cadre de la methode peuvent 
, contenir les informations de variation de 1' amplitude avec le deport ou 
offset : donnees avant sonmiation (« prestack »), donnees liees a la 
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technique connue de rhomme de I'art dite « intercept - gradient » 
extraites de T analyse de cette variation, etc.) ou ne pas les contenir : 
donnees apres sommation (« poststack »). 

La methode s' applique notaniment dans le cadre de T exploration 
5 petroliere pour obtenir des representations quantitatives et des images 
plus fines de la structure ou configuration d'une zone souterraine et elle 
facilite la caracterisation des gisements d'hydrocarbures. 

Etat de la technique 

Des methodes d' inversion ont ete developpees, permettant de relier 
10 a des donnees sismiques, a des donnees mesurees in situ et a des 
observations ou interpretations, une grandeur physique d'un milieu 
heterogene telle que T impedance qu'il presente relativement k des ondes 
elastiques. 

On part d'un modele a priori^ (ID, 2D ou 3D) d'un ou plusieurs 
15 parametres physiques du milieu construit par interpretation des donnees 
connues, et Ton definit des operateurs de covariance decrivant les 
incertitudes sur ce modele a priori et aussi les incertitudes sur les dormees 
connues. On construit des sismogranunes synthetiques qui constituent la 
^ reponse du modele. On les compare aux sismogranunes reels obtenus par 
20 exploration sismique et on minimise de proche en proche les ecarts 
constates entre eux selon des normes associees a des operateurs de 
covariance choisis en fonction d'une evaluation des incertitudes sur le 
modele a priori et les enregistrements sismiques. 
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Les differentes approches connues different notamment par le 
nombre de parametres physiques decrivant le milieu, la dimension du 
modele a priori, et la richesse des possibilites offertes dans la prise en 
compte des parametres d' incertitude. EUes consistent a obtenir le modele 
5 optimal qui minimise une fonction objectif, somme d'une fonction 
objectif dite sismique, et d'une fonction objectif dite geologique. 

La fonction objectif sismique represente Tecart (au sens de la 
norme induite par Toperateur de covariance relatif aux donnees) entre les 
donnees observ6es et les donnees synthetiques. Ces demieres resultent de 
10 la prise en compte d'un modele de propagation plus ou moins complexe. 

La fonction objectif geologique mesure Tecart au module a priori 
(au sens de la norme induite par I'operateur de covariance relatif au 
milieu). 

Parmi ces m^thodes, on pent citer celles d^crites par Brae J. et al; 
15 "Inversion with A Priori Information: an Approach to Integrated 
Stratigraphic Interpretation"; in Sheriff R.E.Ed., Reservoir Geophysics, 
Soc. ExpL Geophys., Investigations in Geophysics 7). 

Une autre m6thode de modelisation en 2D d'une grandeur physique 
par inversion stratigraphique est egalement decrite dans le brevet EP 
20 354.1 12 dudemandeur. 

Par le brevet FR 2.765.692 du demandeur, on connait aussi une 
methode du meme type pour modeliser en 3D un parametre ou grandeur 
physique telle que T impedance d'un milieu heterogene, avec construction 
d'un modele a priori k partir d'un modele geometrique en 3D comprenant 
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plusieurs volumes feuilletes (unites sedimentaires) et aussi de valeurs 
connues de cette grandeur en differents points du milieu, la selection 
d'un modele de covariance le long des feuilles des differents volumes et 
la formation d'un modele optimal par un processus d' inversion a partir 
5 du modele a priori. 

Les methodes precede ntes ne permettent de modeliser que la 
variation d'un seul parametre physique dans le milieu considere. De plus, 
les parametres que Ton utilise pour decrire les incertitudes sur le modele 
a priori (parametres definissant I'operateur de covariance) sont constants 
10 dans tout le domaine ou par unite geologique. 

Des methodes existent qui perniettent de modeliser plusieurs 
grandeurs ou parametres physiques dans un nndlieu. On en trouve decrites 
notamment par : 

Simmons et Backus, 1996 « Waveform-based inversion and AVO 
15 prediction-error », in Geophysics, 61, p. 1575-1588 ; ou par 

Pan et al., 1994 « An Integrated Target Oriented Prestack Elastic 
Waveform Inversion : Sensitivity, Calibration, and Application », in 
Geophysics, 59, 9, p. 1392-1404. 

Mais les modeles a priori consideres dans le cadre de ces methodes 
20 connues, sont de type ID (ils ne dependent que de la profondeur), et 
chaque groupe de traces relatif a une position laterale donnee du milieu 
(CDP), est traite independamment d'une position laterale k une autre. 
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La methode selon rinvention 

La methode selon I'invention permet d'etendre Tapproche 

developpee dans les deux brevets precedemment cites, au cas de 

plusieurs grandeurs ou parametres physiques en offrant de plus la 

5 possibilite de faire varier en tout point du milieu les parametres 
dUncertitude decrivant les incertitudes sur le modele a priori, moyennant 
certaines conditions. 

La m6thode a pour objet de determiner un modele optimal en 2 ou 3 
dimensions (2D ou 3D), representatif de la variation dans un milieu ' 
10 h6t6rogene de plusieurs parametres physiques a partir de donnees 
. connues obtenues par differents moyens : donnees enregistrees 
correspondant a des ondes renvoyees par les discontinuites du milieu en 
reponse ^ des ondes qui s'y propagent et a partir d'un modele a priori (2D 
ou 3D) de ce milieu heterogene, tout en prenant en compte les 
15 incertitudes sur les donnees enregistrees ainsi que celles sur le modele a 
priori. 

La methode selon Tinvention permet d'obtenir un modele optimise 
representatif en 2D ou 3D d'un milieu heterogene tel que le sous-sol, 
decrit par au moins un ou plusieurs parametre(s) ou grandeur(s) 
20 physique(s), i partir de donnees enregistrees correspondant a des ondes 
renvoy6es par le milieu et capt6es par des r^cepteurs couples avec le 
milieu, en reponse a des ondes qui y sont transmises par une source, et 
d' information a priori incluant des donnees mesurees in situ. 
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EUe comporte la construction d'un modele geometrique en 2D ou 3D 
decrivant des lignes ou surfaces de correlation, la construction d'un 
modele a priori en 2D ou 3D decrit par un ou des parametre(s) 
physique(s), a partir de ce modele geometrique et des donnees mesurees 
5 in situ de ces parametres physiques faites a differents points du milieu 
heterog^ne, et la formation d'un modele optimal par inversion des 
donnees enregistrees en prenant en compte toutes les informations a 
priori. 

EUe est caracterisee en ce que Ton choisit un operateur de modelisation 
lO des incertitudes sur le modele a priori, de telle sorte qu'il autorise la prise 
en compte d'un degre de confiance sur le modele a priori variable en tout 
point. 

On choisit par exemple uii modele de covariance exponentiel qui 
peut etre variable selon la localisation dans le milieu. Le modele de 
15 covariance selectionne peut etre selon les cas a ID, a ID le long des 
lignes de correlation, a 2D, a 2D le long des surfaces de correlation ou 
bien encore a 3D. 

Suivant un ou plusieurs modes de mise en oeuvre, le modele de 
covariance exponentiel est relatif a plusieurs parametres. 

20 Le modele etant decrit par plusieurs parametres physiques, on 

modifie par exemple T operateur de covariance de fa9on a definir des 
ecarts par rapport au modele a priori pour d'autres parametres physiques, 
de telle sorte que ces 6carts soient independants en chaque point donne 
du milieu. 
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Le modele de covariance est par exemple ariisotrope, avec des 
longueurs de correlation diff^rentes selon diff^rentes directions dans le 
milieu. 

Suivant un mode de mise en oeuvre, la formation du modele optimal 

5 a plusieurs parametres physiques est obtenue en minimisant une fonction 
objectif globale comprenant un terme relatif aux donnees enregistrees, 

proportionnel au carre de la norme de Tecart entre des donnees 
synthetiques et les donnees enregistrees, et d'un terme relatif au milieu, 
constitue dans le cas d'un modele a 2D respectivement 3D) d'une part du 

10 carre de la norme de Tecart par rapport au modele a priori exprime par 
les dits autres parametres physiques, et d' autre part, du carre de la norme 

L^du gradient de cet ecart, projete de ce gradient sur la tangente a la 
ligne de correlation locale; (respectivement sur le plan tangent a la surface 
locale de correlation). 

15 Dans le cas d'un modele a priori en 3D, on peut former le terme 

relatif au milieu a partir du carr6 de la norme de T ecart au modele a 
priori mesur6 relativement aux dits autres parametres physiques, de la 

norme du gradient de cet ecart, et du carre du Laplacien de T ecart. 

La methode selon I'invention peut etre appliquee a la modelisation 
20 d'une zone du sous-sol, les mesures in situ etant obtenues a differentes 
profondeurs dans un ou plusieurs puits au travers de la zone, les ondes 
emises etant des ondes elastiques. 
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La methode selon T invention peut etxe aussi appliquee a la 
modelisation d'une zone d'un nailieu ou les ondes emises sont des ondes 
electromagnetiques , 

L'interSt de I'approche selon r invention est : 

5 - d'une part d'estimer plusieurs parametres physiques decrivant le 

milieu heterogene, ce qui a pour consequence de mieux caracteriser ce 
milieu, et 

- d' autre part de prendre en compte des informations a priori tres 
riches, tant au niveau de la dimension du modele a priori que des 

10 incertitudes le concemant. 

Presentation des figures 

D'autres caract6ristiques et avantages de la methode selon 
r invention apparaitront a la lecture de la description ci-apres et de 
resultats bbtenus sur un exemple de mise en oeuvre, en se referant aux 
15 dessins annexes ou : 

- les Fig. 1 et 2 montrent les resultats obtenus en impedance P et 
impedance S par une premiere inversion a partir des seules donnees,PP, 
c'est a dire sans informations a priori ; 

- les Fig. 3 et 4 representent les impedances P et S de ce reservoir 
20 realiste apres filtrage dans la bande de frequence de la sismique ; 

- les Fig. 5 et 6 montrent les imp6dances P ef S du module a priori 
calcule par interpolation des donnees de puits le long des lignes de 
correlation ; et 
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- les Fig. 7 et 8 representent les resultats en Ip et Is pour une 
inversion avec information a priori ou la longueur de correlation dans 
Tunite Uj est fixee a 500m dans T unite Ul, et a 50m dans les unites U2, 
U3,etU4; 

5 Description detaillee 

On va considerer a titre d'exemple une zone souterraine que Ton 
veut imager en faisant apparaitre comment varient dans ceue zone 
plusieurs parametres physiques tels que les impedances relativement aux 
ondes P et S, la masse specifique etc. On dispose d' informations a priori 

10 permettant de former un modele a priori a 2D ou 3D de la zone decrite 
par ces parametres physiques. Ces informations a priori sont par exemple 
derivees d'enregistrements sismiques obtenus classiquement en captant 
au moyen de recepteurs sismiques Ri couples avec la formation 
(prospection sismique terrestre) ou remorques en immersion par un 

15 navire (prospection sismique marine), les ondes renvoyees par les 
discontinuites en reponse a des ondes sismiques emises par une source S. 
Les informations a priori sont aussi obtenues en realisant des diagraphies 
dans un ou plusieurs puits W fores au travers de la formation. On dispose 
aussi generalement d' informations structurelles d'ordre g^ologique sur la 

20 formation. 

On realise T optimisation de ce modele a priori par inversion suivant 
le schema general dej^l utilise dans les brevets pr6cit^s consistant ^ 
minimiscr les ecarts entre des enregistrements synthetiques (constituant 
la reponse du modele en cours d' optimisation) et les enregistrements 
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reels dont on dispose. A cet effet, on minimise une fonction-objectif 
globale a deux termes : un teraie sismique et un tenne geologique, au 
moyen d'une m6thode de gradient conjugue par exemple. 

Le premier terme de la fonction-objectif, qui est le terme sismique, 
5 mesure Tecart entre les donnees observees et les donnees synthetiques au 
sens de la norme associee a T inverse de Toperateur de covariance, note 
Q, qui decrit les incertitudes sur les donnees. Le second mesure Tecart 
entre le modele courant et le modele a priori, au sens de la norme 
associee Tinverse de I'operateur de covariance, note c„ , qui decrit les 
10 incertitudes sur les modeles. En minimisant cette fonction-objectif 
globale, on peut obtenir en une seule passe un modele multiparametre, 
compatible k la fois avec les donnees sismiques observees et le modele a 
priori multiparametre. 

La methode selon T invention permet T application de ce schema 
15 general de modelisation a une formation decrite par plusieurs parametres 
physiques. A cette fin, on va utiliser une fonction-objectif globale dont 
les deux termes sismique et geologique vont etre definis plus en detail ci- 
apres. 

La fonction objectif sismique 

20 La fonction objectif sismique mesure Tecart entre chaque trace 

sismique et la trace sismique qui lui est liee, au sens de la norme associee 
a r inverse de Toperateur de covariance On suppose que le bruit 
sismique est decorrele d'echantillons en echantillons. En consequence, 
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on choisit un operateur de co variance diagonal decrivant T incertitude sur 
les donnees, dont la fonction de covariance est : 

ou a] est la variance de T incertitude sur les donnees, c'est-a-dire la 
5 variance du bruit sismique, et bd{xj,m) represente Tecart entre les 
donnees synthetiques et les donnees observees au point de coordonnee 
{x,t,w) (r indiquant le temps, x la position laterale, et m Toffset (c'est a 
dire la distance source-recepteur)). 

Cet operateur est inversible, et la norme associee a son inverse 
10 permet de definir Texpression de la fonction objectif sismique par : 

\\d-d"%~jjj{d-d'''''y{x,t,XD)dxdtduy 

Dans cette 6quation, d*'^' represente les donnees observees, et les 
donn6es synthetiques d sont calcul6es h partir de la donn6e d'un modele 
de sous-sol m et d'un operateur de modelisation G, selon I'expression 
15 suivante : 

d-^G{m) 

La fonction objectif geologique 

On sait que Tinversion 2D implique un modele de covariance ID le 
long des lignes de correlations, et T inversion 3D, un modele de 
20 covariance 2D le long des surfaces de correlations. Pour simplifier 
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Texpose qui suit, on se place dans le cas a 2 dimensions, sachant que les 
equations suivantes se generalisent facilement au cas 3D, 

Soit 5m(p)=m(p)-mP^'°''(p) la fonction ecart au modele a priori au 
point p. Conune fonction de covariance de Toperateur Cm, on choisit la 
fonction : 

cov(5m(p).5m(pO) = Q(p)D''=(p)D'-ip-X»(pOexp(-||A(p)-X(p'|) 

ou A est une fonction de derivee partout non nulle,.et Q et D deux 
fonctions matricielles, telles que pour tout point p, Q(p) et D(p) soient 
deux matrices respectivement orthonormale et diagonale. Ces deux 
matrices proviennent de la decomposition en vecteur propre/valeur 
propre de la matrice suivante: 
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ou k est le nombre de parametres elastiques decrivant le milieu. 

Cette matrice est la matrice de variance/covariance decrivant au 
pomt p, r incertitude sur les parametres (par la variance ai ), et les 
correlations de ces incertitudes (par le coefficient de correlation Py). 

En introduisant les nouveaux parametres definis par 
5m = D''''P'5m 
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et en effectuant le changemenf de coordonnees p = A{p), la 
fonction de covariance precedente devient : 

cov(6m(^),5m(p'))= exp(--jlp- p|) 

II est a la portee de Thomme de Tart d'inverser un operateur ayant 
une telle fonction de covariance. 

Le terme geologique que I'on utilise ici est la norme de Sm associee 
a r inverse de T operateur de covariance que Ton vient de decrire. 

Dans le cas ou Ton desire que les incertitudes soient decrites le long 
des lignes de correlation en 2D, alors en notant s Qt x les coordonnees 
curvilignes le long des lignes de correlations, la direction de s etant 
tangente au lignes de correlation. 

cette demarche se generalisant aisement au cas 3D 

15 Dans le cas par example, ou Ton peut supposer que A(x) = x/X, X 

6tant la longueur de correlation des incertiuides le long des lignes de 
corrdlation, le terme geologique de la m6thode est : 

ou s[x) est le vecteur tangent a la ligne de correlation passant au 
20 point X, et 5m est definit par : 

5m = D'^'^P'6m 


( 


5 
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Ainsi, le terme geologique que Ton a choisi permet de prendre en 
compte des parametres d' incertitudes (variances, longueurs de 
correlations, et coefficient de correlation des incertitudes inter- 
parametres) qui peuvent varier en tout point du domaine. 

5 Validation de la methode 

II s'agit d'un exemple synthetique deriyee d'un cas reel etudie suite 
a des operations de prospection sismique marine, oii la zone reservoir 
travefsee par deux puits Wl, W2 est incluse dans un milieu encaissant 
ID. 

10 Des donnees sismiques synthetiques ont ete simulees, auxquelles on 

a ajoute un bruit aleatoire tel que le rapport bruit sur signal soit egal a 
80%. Les ondes PP reflechies sur le reservoir constituaient les donnees 
dites « observees ». On dispose de 240 enregistrements en couverture 
multiple (CDP) avec 50 valeurs de deport (« offset ») par couverture 

15 multiple avec un deport minimum de 50m, et un deport maximum de 
1275m. 

L' information geomeuique du modele a priori a ete construite a 
partir de 3 horizons interpretes sur le resultat en impedance PP, sans 
informations a priori. Ces trois horizons definissaient quatre unites 
20 geologiques Ui, U2, U3, et U4. Des lignes de correlations ont ete definies 
en tout point d'apres une hypothese stratigraphique de conformite suivant 
ces trois horizons limites. 

L' information en impedance P et S aux deux puits Wl, W2 aux 
abscisses respectives x=250m et x=5750m, a ensuite ete introduite. 


On a choisi une matrice de covariance exponentielle le long des 
lignes de correlation, avec une longueur de correlation X. La variance de 
rincertitude sur Tecart au modele a priori etait prise egale a 10% pour 
r impedance P, et a 20% pour V impedance S. 

Ces resultats avec information a priori compares avec ceux obtenus 
sans information a priori illustrent bien, d'une part en quoi T information 
a priori permet d'ameliorer 1' image obtenue en impedance S et, d' autre 
part en quoi T information a priori permet de diminuer les incoherences 
entre les traces, visibles sur les Figures 1 et 2. 
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REVENDICATJONS 

1) Methode pour obtenir un modele optimise representatif en 2D 
ou 3D d'un milieu heterogene tel que le sous-sol, decrit par au moins un 
ou plusieurs parametre(s) ou grandeur(s) physique(s), a partir de donnees 

5 enregistrees correspondant a des ondes renvoyees par le milieu et captees 
par des recepteurs couples avec le milieu, en reponse a des ondes qui y 
sont transmises par une source, et d' information a priori incluant des 
donnees mesurees in situ, la methode comportant la construction d'un 
modele geometrique en 2D ou 3D decrivant des lignes ou surfaces de 

10 correlation, la construction d'un modele a priori en 2D ou 3D decrit par 
un ou des parametre(s) physique(s), a partir de ce modele geometrique et 
des donnees mesurees in situ de ces parametres physiques faites a 
diff^rents points du milieu heterogene, et la formation d'un modele 
optimal par inversion des donnees enregistrees en prenant en compte 

15 toutes les informations a priori, caracterisee eri ce que : 

- on choisit un operateur de modelisation des incertitudes sur le modele a 
priori, de telle sorte qu*il autorise la prise en compte d'un degre de 
confiance sur le modele a priori variable en tout point. 

2) Methode selon revendication 1, caracterisee en ce que Ton 
20 selectionne un modele de covariance exponentiel 

3) Methode selon revendication 2, caracterisee en ce que Ton 
choisit un modele de covariance exponentiel variable selon la localisation 
dans le miheu. 
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4) Methode selon Tune des revendications precedentes, caracteris6e 
en ce que le modele de covariance selectionne est a ID ou a ID le long 
des lignes de correlation, ou a 2D, ou a 2D le long des surfaces de 
correlation ou encore a 3D. 

5 5) Methode selon Tune des revendications precedentes, caracterisee 

en ce que le modele de covariance exponentiel est relatif a plusieurs 
parametres. 

6) Methode selon Tune des revendications precedentes, caracterisee 
en ce que, le modele etant decrit par plusieurs parametres physiques, on 

10 modifie I'operateur de covariance de fa9on a definir des ecarts par 
rapport au modele a priori pour d'autres parametres physiques, de telle 
sorte que ces ecarts soient independants en chaque point donne du milieu. 

7) Methode selon Tune des revendications precedentes, caracterisee 
en ce que le modele de covariance est anisotrope, avec des longueurs de 

15 correlation differentes selon differentes directions dans le milieu 

8) Methode selon Tune des revendications precedentes, caracterisee 
en ce que la formation du modele optimal est obtenue en minimisant une 
fonction objectif globale comprenant un terme relatif aux donnees 

enregistrees, proportionnel au carre de la norme de Tecart entre des 
20 donnees synthetiques et les donnees enregistrees, et d'un terme relatif au 
milieu, constitue dans le cas d'un modele a 2D (respectivement 3D) d'une 

part du carre de la norme de I'ecart par rapport au module a priori, et 
d'autre part, du carre de la norme L^du gradient de cet ecart, projete sur 
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la tangente a la ligne de correlation locale (respectivement sur le plan 
tangent a la surface locale de correlation). 

9) Methode selon rune des revendications precedentes, caracterisee 
en ce que la formation du modele ojptinial a plusieurs parametres 

5 physiques est obtenue en minimisant une fonction objectif globale 
comprenant un terme relatif aux donnees enregistrees, proportionnel au 

carre de la norme de I'ecart entre des donnees synthetiques et les 
donnees enregistrees, et d'un terme relatif au milieu, constitue dans le cas 
d*un modele a 2D (respectivement 3D) d'une part du carre de la norme 

10 l2 de Tecart par rapport au modele a priori exprime par les dits autres 

parametres physiques, et d'autre part, du carre de la norme L^du gradient 
de cet ecart, projete sur la tangente a la ligne de correlation locale 
(respectivement sur le plan tangent a la surface locale de correlation). 

10) Methode selon Tune des revendications - precedentes, 
15 caracterisee en ce que dans le cas d'un modele a priori en 3D, le terme 

- relatif au milieu est proportionnel au carre de la norme de Tecart au 

modele a priori, de la norme du gradient de cet ecart, et du carr6 du 
Laplacien de T ecart. 

11) Methode selon Tune des revendications precedentes, 
20 caracterisee en ce que dans le cas d'un modele a priori en 3D, le terme 

..relatif au milieu est proportionnel au carre de la norme de Tecart au 
modele a priori mesure relativement aux dits autres parametres 
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physiques, de la norme du gradient de cet ecart, et du carre du 
Laplacien de Tecart. 

12) Methode scion I'une des revendications precedentes, 
caract^risee en ce que le milieu heterogene est une zone du sous-sol, les 

5 mesures in situ sont obtenues a differences profondeurs dans un ou 
plusieurs puits au travers de la zone. 

13) Methode selon Tune des revendications precedentes, 
caracterisee en ce que les ondes emises dans le milieu sont des ondeis 
elastiques. 

10 14) Methode selon Tune des revendications precedentes, 

caracterisee en ce que les ondes ennises dans le niilieu sont des ondes 
electromagnetiques. 
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